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We confront the new physics models with extra iso-singlet down quarks in the new CP violation 
experimental era with sin (2/3) and e'/e measurements, K'^ tv^uu events, and Xs limits. The 
closeness of the new experimental results to the standard model theory requires us to include full 
SM amplitudes in the analysis. In models allowing mixing to a new isosinglet down quark, as in 
Ee, flavor changing neutral currents are induced that allow a Z" mediated contribution to B — 13 
mixing and which bring in new phases. In {p,rj), (a;s, sin (7)), and (a;^, sin (2(/>s)) plots we still 
find much larger regions in the four down quark model than in the SM, reaching down to 77 « 0, 
< sin (7) < 1, —.75 < sin (2a) < 0.15, and sin (20s) down to zero, all at la. We elucidate the 
nature of the cancellation in an order four down quark mixing matrix element which satisfies the 
experiments and reduces the number of independent angles and phases. We also evaluate tests of 
unitarity for the 3x3 CKM submatrix. 

PACS numbers; 11.30.Er,12.15.Hh,12.15.Mm,12.60.-i,14.40.Nd 
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I. INTRODUCTION 



The "new physics" class of models we use are those 
with extra iso-singlet down quarks, where we take only 
one new down quark as mixing significantly. An example 
is Eg, where there are two down quarks for each gener- 
ation with only one up quark, and of which we assume 
only one new iso-singlet down quark mixes strongly. This 
model has shown large possible effects in B — B mixing 
phases [Q. The new B factory results on sin (2/3) in the 
SM range, the e'/e experimental convergence, the new 
if+ TT+i/i/ result, the Ams limits near the SM predic- 
tion, and other new measurements require a finer analysis 
and a potential challenge to new physics models. In this 
paper we include the full SM contributions as well as 
the new physics contributions from the iso-singlet down 
quark model to jointly analyze the constraints from all of 
these experiments, as well as other flavor changing neu- 
tral current (FCNC) limits and SM CKM matrix element 
constraints. 

In models allowing mixing to a new iso-singlet down 
quark (as in Eg) flavor changing neutral currents are in- 
duced that allow a mediated contribution to B — B 
mixing and which bring in new phases ^, ^ . In {p, rj) , 
(xs, sin (7)), and (xs, sin (2(/)s)) plots we still find much 
larger regions than in the SM, reaching down to 77 « 0, 
< sin (7) < 1, and sin(2(/)s) down to zero (below the 
SM range), all at la limits. The nature of the cancella- 
tion in a fourth down quark matrix element Vm to satisfy 
the experiments is elucidated. We also establish ranges 
for the new mixing elements to the new iso-singlet down 
quark, and make a simple estimate of the lower mass 
limit of the new down quark. 
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In Section II we introduce the scenario with more down 
quarks as in Eg, truncate it to one extra down quark, in- 
troduce the 4x4 mixing matrix, and apply it to -B — i? 
mixing. Section III presents the CP violating Bd and 
Bs decay asymmetries, and Bs mixing, including the 
FCNC tree diagram additions. Section IV presents the 
full SM contributions as well as the four down quark 
model (FDQM) amplitudes for the CP violating and 
FCNC K meson experiments that are used. Section V 
presents the joint chi-squared analysis and results for the 
SM and FDQM model for the various plots listed above. 
Section VI presents the sizes or limits on the matrix ele- 
ments, mixing angles, phases, FCNC couplings and uni- 
tarity quadrangles. Section VII lists the conclusions and 
projects what the next down quark mass limit might be. 



II. ISO-SINGLET DOWN QUARK MIXING 
MODEL 



Groups such as Eg with extra SU(2)l singlet down 
quarks Q give rise to fiavor changing neutral currents 
(TCNCVthrough the mixing of four or more down quarks 
lUi 01 H • The initial quarks of definite weak isospin in 
Eg, for each generation are: the left handed iso-doublet 
(M°^,d°^), their right handed iso-singlets and 
and the yet to be found iso-singlet pairs D^j^ and 1?°^. 

We can take the initial up quark matrix to be the mass 
eigenstates, so m° = Ui, giving = /3X3. The down 
quarks {d^,D^), which correspond to the same genera- 
tions as Ui, mix to form mass eigenstates (c?i, Di) via the 
matrix ^''(6 x 6), where d"^ = Vj^djL- The weak inter- 
action charged current matrix is then U ~ Y"^ x y*, the 
3x6 matrix that is the upper three rows of V^. The 
lower three rows of are the three linear combinations 
of {di, Di) that are the iso-singlet which cannot couple 
to up quarks by the weak interactions. 
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We truncate the V matrix to the 4x4 matrix using 
only the D quark that mixes most (and dropping the su- 
perscript d on V^), giving the Four Down Quark Model 
(FDQM). Calling the new down quark mixture D, the 
weak charged currents of Z? to it, c, and t quarks are 
VtD = S34, VcD = 5246-**^*, and VuD = si4e~'^'\ which 



are in the fourth column. The fourth row gives the linear 
combination that is the initial iso-singlet D^. The com- 
plete 4x4 mixing matrix was given previously |9[ p^ . 
The leading terms in the 4x4 down quark mixing matrix 
with 6 angles and 3 phases are 




d s b 

u / C12C34 S12C34 sise"*^'^ 

T r _ C I — S12 1 S23 ■^:^4'- I 

H (si2S23"Sl3e^'") -S23 1 .^.1' 

4 V V^d 14. V^h 

where, to leading order in new angles, 

Vld = -si4e-'^"+S24e-'^^^si2-S34(si2S23-si3e-'^i^), (2) 

Yl = _s24e-'524_5^^g-»5i4^^2 + g3^(^23 + 5125136"'^^^), (3) 

= -S34-S24e-**^^S23-si4e-'^"si3e^*^^ (4) 



A. FCNC in Z° Couplings From Extra Iso-singlet 
Down Quarks 

The FCNC amplitudes are given in terms of the mix- 
ings p4i to form the iso-singlet down quark by [pi 



t^y = ^l^^3 for « 7^ i- 



(5) 



The FCNC couplings of the down quarks to the are 
then given by 



FCNC 



2 sin 9w cos 6w 



Uijd.,Ll^djLZ^. (6) 



The flavor changing neutral currents are |7[ [sj —Usd — 
V^Vm , -Usb = V^Vib, and ^Um = Vl^VAd- 

The diagonal neutral current couplings are reduced 
in strength by the amplitudes into the iso-singlet down 
quarks, becoming 



'"NC — 



2 sin6'vi/ cos Bw 



The FCNC with tree level mediated exchange may 
contribute part of B\ — B'^ mixing and of B'^ — B^ mixing, 
and the constraints leave a range of values for the fourth 
quark's mixing parameters. As shown in Fig. ^, B^ — B^ 
mixing may occur by the b—d quarks in a B^ annihilating 
to a virtual Z through a FCNC with amplitude Um , and 
the virtual Z then creating b — d quarks through another 
FCNC, again with amplitude Utd, which then becomes a 
Bd meson. 

If the FCNC amplitudes are a large contributor to the 
Bd — Bd mixing, they introduce three new mixing angles 
and two new phases over the standard model (SM) into 
the CP violating B decay asymmetries. The size of the 




FIG. 1: The SM second order weak box diagram plus the 
double FCNC vertex tree diagram with an intermediate Z'^ 
for Bd - Bd mixing 



contribution of the FCNC amplitude Udb as one side of 
the unitarity quadrangle is less than 0.15 of the unit base 
l^cdKbl at the l-CT level (see Section VI), but we have 
found [1,1,0, 1 that it can contribute as large an amount 
to Bd — Bd mixing as does the standard model. The new 
phases can appear in this mixing and give total phases 
different from that of the standard model in CP violating 
B decay asymmetries @, ||, ||, |l2| . 

For Bd — Bd mixing with the four down quark induced 
b—d coupling, Udb, we have [0[ 



Xd 



{2GF/iV2)BBflmBriBTB \U, 



2 

std—db 



u, 



db\ 



(8) 



where with yt — ml jrr? 



C/2,_,, = (a/(4^sin2 0w))2/t/2(^/t)(V;rfl46)^ (9) 

and Xd = I^rriB^I^ Bd — TBd^'^Bd- In order to compare 
magnitudes, in the SM, C/^td-dh = 0.50x 10"^(l-p+i77)^. 

The CP violating decay asymmetries depend on the 
combined phases of the B^ — B'^ mixing and the b quark 
decay amplitudes into final states of definite CP. Since 
we have found that Z mediated FCNC processes may 
contribute significantly to B^ — B^ mixing, the phases 
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of Udb would be important. The FCNC amplitude Udb 
to leading order in the new angles is 

Udb = (-S34-S24S23e**^^) 

[s^^V^d + si4e-^''" - S24e-*^^^si2), (10) 
where Vtd ~ (S12S23 - siae*'^"), and Kb = siae"**^". 



III. MIXING AND CP VIOLATING DECAY 
ASYMMETRIES IN THE FOUR DOWN QUARK 
MODEL 

With new additive contributions to CP violating decay 
asymmetries, the asymmetries are no longer sines of SM 
unitarity triangle angles. However, they are still sines of 
the overall phases of the amplitudes in the asymmetries. 
We analyze the FDQM with the present data, and also 



show projected results for three different sin(2Q;) values of 
— 1.0, 0, and +1.0, which are allowed under the FDQM, 
although sin(2a) — ±1 are not allowed by the SM at 2a. 



A. sin (2/3) and sin (2a) 

In the four down quark model we use "sin (2q!)" and 
"sin (2/3)" to denote results of the appropriate Bd de- 
cay CP violating asymmetries, but since the mixing am- 
plitudes are superpositions, the experimental results for 
these asymmetries are not directly related to angles in a 
triangle. Being imaginary parts of pure complex expo- 
nentials, they are sines of phase angles. The asymmetries 
with FCNC contributions included are (for B mixing to 
B before decay) 



sin (2/3) = Abo^ 



Im 



{Uld- 



db + uj,) {v:,v,s 



) {v:,Vud)' 



db 



ulb\ 



(11) 



sin {2a) = -A^ 



-Im 



' {uld-db + ulb) (KTbKd) " 
.\uld-db + ul\ {VubV:,)_ 



(12) 



with U^fd-db defined in Eq. (^). The same mixing phase 
occurs in both asymmetries, times the squares of the dif- 
ferent decay phases. We take the Moriond 2002 results 
for sin (2/3) from Babar ^ and Belle ||l|, || to give the 
weighted average sin(2/3) = 0.78 ± 0.08. 



the Z-mediated annihilations and exchange amplitudes 
as well. Here we avoid the full theoretical uncertainty on 
Bb/b, by taking the ratio of to Xd, which is better cal- 
culated theoretically, and in which we have also included 
the FCNC with Z° exchange 



B. sin (7) 



Am, 



1.35 Xd 



\^std-sb 



"-^std-db 



u!b\ 



(15) 



In the four down quark model, what we mean by 
"sin (7)" is the result of the experiments which would 
give this variable in the SM (l|, |l|l, as in B° D+R-'. 
Here, the four down quark model involves more compli- 
cated amplitudes, and "sin (7)" is not simply sin((5i3): 



sin(7) = Im 



where 



(uL-sb + ul) (v;:,K.) (k;k.) 



\^'std-sb 



cb' 



, (13) 



U. 



std—sb 



{a/{A7rsmew'))ytf2iyt){V*Vtbf. (14) 



In the SM, C/2 



10 X 10-*^. 



C. The "Frequency" of Bs Oscillations, Xs 

In the four down quark model, x^ is no longer the sim- 
ple ratio of two CKM matrix elements, but now involves 



We now include the amplitude method analysis of LEP 
with SLD to assign a Ax^ for each Am^ calculated in the 
angular parameter grid |18|. 



D. The Bs Decay Asymmetry, sin(2</)s) 

In the standard model, Bs mixing involves {Vf*Vtb)'^ 
which is almost exactly real, and the leading decay pro- 
cess of 6 — > CCS has no significant phase from the decay 
which is proportional to Vj,. Thus almost no CP vio- 
lating phase develops in the most likely Bg decays. This 
occurs in the decays Bs J/'^<f>, Bs — > DfDj, and 
Bs J/^Ks- The near vanishing of this asymmetry 
is a test of the SM |^]. Below, we will find a strange 
twist on this, since the FDQM will include a range that 
includes values smaller than the SM range, and does not 
exceed it. In the SM the angle (j)s is the small angle in the 
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b — s unitarity triangle, and its non-zero value indicates 
CP violation. 

In the four down quark model, the CP violating Bs 
decay asymmetry is (for the mixing to J/iJj cj) or 
without the final Ks) 



sin [2(f) s) 



-Im 



\^std~sb 



(16) 



which includes the double FCNC exchange propor- 
tional to U^g. Because of the additional flavor changing 
term, in the four down quark model, the angle given by 
the above asymmetry will not generally be an angle in a 
triangle. 



IV. FOUR DOWN QUARK MODEL 
AMPLITUDES IN KAON EXPERIMENTS 

A. FCNC as an addition to Penguin plus Box 
Amplitudes 

Since CP violation and FCNC experiments with K 
mesons are approaching the SM range and also limit 
FCNC amplitudes, we now include the full SM ampli- 
tudes with the FCNC Z° exchange amplitudes as well. 
The K meson experiments are e, -n^vv, y^^ij 

and we now add the recent and fairly well determined re- 
sults for Re(e/e'). 

We use the amplitudes determined by Buras pO| . 
In order to reconcile the notation between us and Buras 
and Silvestrini pl| |, we relate their Zds to our Ugd by 
taking 



Zds — — 



V2C 



sdi 



(17) 



as implied by their definitions in Lagrangians. 

In the following formulas, Bq is the AS" — 1 box am- 
plitude. So is the AS* — 2 box amplitude, Cq is the Z'^ 
Penguin amplitude, Dq is the off-shell photon penguin 
(with Dq being the on-shell amplitude), and Eq is the 
off-shell gluon penguin {Eq being on-shell). Gauge inde- 
pendent combinations are 



^0 = 
Yo = 

Zo = 



Co - 4Bo 
Co — Bo 

Co + \do 



(18) 
(19) 

(20) 



For mt — 170 GeV and rric — 1.25 GeV, for exam- 
ple, these quantities are So{xt) — 2.46, Sq{xc) — Xc, 
So{xc,Xt) = 0.0022, Xo = 1.57, yo = 1-02, Zq = 0.71, 
Eq 0.26, D'q = 0.38, and E'q = 0.19. 

The FCNC Z^ exchange with amplitude Uds can be 
added to the d~ s Penguin amplitude with the Z" by the 
substitution 



AtCo [xt) -> XtCo {xt) 



to obtain the SM plus FCNC result. 



Usd, 



(21) 



B. Indirect CP Violation in Epsilon 

In K — K mixing, the small indirect CP violation is 
given through |e| |l9] 



1 



|ImMi2|, 



V2AMk 

where we include the substitution from Eq. Ell 



(22) 



C^ - 

Mi2 = —^Fl,BKmKoM^[XfTjiSo{xc) + XfmSo(xt) + 2X*\lf]3So{xc,xt)] 

-^Fl,Bj,m^oU!,. (23) 



The short distance QCD corrections factors in NLO are 
@ ryi = 1.38 ± 0.20, 772 = 0.57 ± 0.01, and 773 = 0.47 ± 
0.04, and we use Bk = 0.85 ± 0.13. 



C. Direct CP Violation in Re(e7e) 

The direct CP violation in decays, Re(e'/e), has 
received more accurate measurements that are definitely 
non-zero. The average of KTeV [|2| and NA48 ||2^ gives 
Re(e'/e) = 17.3±2.4, where the error has been increased 



by VWJW)- The sum of the SM |T9[ g plus FCNC 
amplitude from Eq. ^ is 

^2 



Re(e7e) = ImAtF^' 



where 



-^=^Imt/.4Px+P.+P.], 

(24) 



-F,' = Po+PxXo{xt) + PYYo{xt)+PzZo{xt)+PEEo{xt). 

(25) 

The P's are functions ofS^^ = 1.0±0.3, Bg^^ 0.8±0.2, 
and = 340 ± 50 MeV ^ . 
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D. K+ ■K+iy-p 

The recent detection of two events in 
produced the experimental resuh [E4[ 



BR, 



expt 



TT^vv has 



(26) 



compared to the SM range || of (0.72±.21) x IQ-i". The 
Poisson probabiUty for the angle parameters is converted 
to a chi-squared form [ p6| which is convolved into the 
total formula. For this experiment using a logarithmic 
prior, with 2 x Uobs = 4 degrees of freedom |26 , the 
addition to is 



= 2 < 71 >= 2 X (2 events) x 



BR, 



calc 



BR, 



(27) 



'CXpt 



The sum of the SM plus FCNC contributions is ob- 
tained by using Eq. pi]) 



BR 



■calc 



(A- 



TT+iyV) = rKBR{K+ -K^e+v) (28) 

^,2 



|K.p27r2 sin^ Ow 



KXl^ + \tX{xt) 



-Usd 



w 



Here, X{xt) = r]^Xo{xt), etc., and [g^ ry^; = 0.994. With- 
out the SM, the contribution of the Z° exchange alone 
with amplitude Usd is — 1.61 x lO^lt/sdP- 



E. Kl 



The short distance weak FCNC contribution to Kl 
fj-IJ, constrains Ke{Uds) and is given from after includ- 
ing Eq. P 



BR {Kl fi'^n') = BR{K 
iKspTT^ sin''6'vy 
X I ReAclWi -^ReAtr(a;t) 



(29) 



\/2GfM2, 



■Ref7, 



sd 



JOINT CHI-SQUARED ANALYSIS OF THE 
SM AND THE FDQM EXPERIMENTS 



FCNC experiments put limits on the new mixing an- 
gles and constrain the possibility of new physics con- 
tributing to B^ — B'^ and B^ — B'^ mixing. Here we jointly 
analyze all constraints on the 4x4 mixing matrix ob- 
tained by assuming only one of the SU(2)l singlet down 
quarks mixes appreciably We use the seven experi- 
ments for the 3x3 CKM sub-matrix elements 1^ , which 
include: those on the three matrix elements Vus, Vut, Vet 
of the u and c quark rows; |e|; Bd — Bd mixing (xd)', the 
new limits on Amg, or Xs', and the new measurements 
for sin (2/3). For studying FCNC, we include Vud and 
Vcd, the bound on B —^ ^^Xg (which constrains b s), 
the two events in K^ ir^vv \12, |6[ and Rb in 
— > bb |l^ (which directly constrains the V46 mixing 
element). FCNC experiments will bound the three ampli- 
tudes Uds, Usb, and Ubd which contain three new mixing 
angles and three phases. We use the mass of the top 
quark as nit ~ 174 GeV. We also add FCNC constraints 
from Kl — > ^J■|J■, now including the large long distance 
error, and the new and more convergent results for e'/e 
from NA48 [|| and KTeV |2|l. 

Related analyses including both SM and FDQM am- 
plitudes in kaon constraints by Barenboim, Botella and 
Vives ^ preceed this work. We have applied a full 
analysis rather than just 95% CL bounds, and have 
included the new, larger and more exact sin (2/3) results, 
as well as new K~^ — > n'^vv, e'/e results, and new and 
full Xs data. We have also included an analysis of the 
4x4 mixing matrix parameters and found a crucial can- 
cellation in one of the matrix elements. 

We use a method for combining the Bayesian Poisson 
distribution for the average for the two observed events 
in K^ TT^vv ^ with the chi-squared distribution 
from the other experiments. This treats the two events 
with a logarithmic Bayesian prior as four degrees of free- 
dom. This gives a total of ten additional experimental 
degrees of freedom for the FDQM. 

In maximum likelihood correlation plots, we use for 
axes two output quantities which are dependent on the 
mixing matrix angles and phases, such as (p, 77), and for 
each possible bin with given values for these, we search 
through the nine dimensional angular data set of the 4x4 
down quark mixing angles and phases, finding all sets 
which give results in the bin, and then put into that bin 
the minimum among them. To present the results, we 
then draw contours at several x^ in this two dimensional 
plot corresponding to given confidence levels. 



Here, Y{xt) = r]yYo, and Q Tyy = 1.012. The long dis- 
tance contribution has been analyzed ||2^ . We make the 
1(7 limit conservatively as the sum of the Icr experimental 
limit plus the Icr long distance estimate |2^. 

From the above K meson formulas, the error formulas 
were generated using Mathematica. 



Standard Model (sin (2a), sin (2^)) Plot 
Present Constraints 



For the SM we take A = Vus as fixed and then use 
the six experiments on the 3x3 CKM matrix elements 
named above with three parameters to give three degrees 
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SM Preser: 1 ci, 90% CL, 2a 




0.5 - 
0.4 V 
0.3 V 

0.2 ^ 
0.- V 

sin(2a) 

sin(2;3) vs. sin(2a) 

FIG. 2: The (sin (2a), sin (2/3)) plot for the standard model 
with contours at Icr, 90%CL, and 2cr with present data. 



FDQM - Present, la, 90% CL, 2a 



0.4 
0.3 

0.2 ^ 
0.1 ^ 

P 

-1 -0.8 -0,6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

sin(2o;) 

sin(2/S) vs. sin(2a) 

FIG. 3: The (sin(2a), sin(2/3)) plot for the FDQM with con- 
tours at Icr, 90% CL, and 2<j. 

of freedom. In the figures we show the contours with 
confidence levels (CL) at values equivalent to Icr, 90% 
CL (1.64ct), and 2a. The new BaBar H and Belle 
|15| average is sin (2/3) = 0.78 ± 0.08. This gives (3 = 
25.6!i^;'4f . From Fig. | for the SM we see that the sin(2/3) 
range is from 0.63 to 0.96 at 2a, centered around the 
experimental average of 0.78. The SM sin(2a) range at 
2cr is from -0.90 to -fO.57. 



B. Four Down Quark Model (sin (2a), sin (2/3)) 
Plots - Present Limits 

In the FDQM analysis including FCNC experiments, 
there are 17 experimental degrees of freedom, minus 9 
parameters, giving 8 remaining degrees of freedom. In 
contrast to the SM, for the FDQM in Fig. |, almost the 
entire region sin(2/3) > 0.48 is allowed at 2ct, and sin(2/3) 
can be as low as 0.55 at Icr. In the FDQM, all values of 
sin(2a) are allowed. In this case, the larger Icr range for 
sin (2/3) than from the direct experimental measurement 
is an effect of including so many experiments in the joint 
fit. 



SM Preser; 1 cj, 90% CL, 2a 




-O.a -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 



p-7] Plot 

FIG. 4: The (p, r?) plot for the SM with three compari- 
son kaon experiments added, with joint contours at la, 
90%CL, and 2cr. The light lines are for the Kaon experiments 
and are described in the text. 



C. Standard Model with Comparison Experiments 
{p,ri) Plot 

Here we depart from the analysis of the SM experi- 
ments alone to show the effects of the additional three K 
meson experiments, namely K'^ — > ir'^vv, ^J^^i, and 

e'/e. While they are not needed in the SM analysis, they 
are included in the FDQM analysis. In this case, there 
are 13 experimental degrees of freedom, minus 4 param- 
eters, giving a net 9 degrees of freedom. Contours are at 
Icr, 90% CL, and 2cr. The three new K meson experimen- 
tal contours for the SM are shown in Fig. ^. For e'/e, the 
lower Icr contour is the horizontal dot-dashed line, where 
the central contour would be a horizontal line at 77 = 1.1. 
For Kl /i/i the solid vertical line sX p — —0.66 is the 
lower Icr contour with the central contour being a ver- 
tical line at p = 0.9, which is not shown. This includes 
conservatively a large and uncertain long distance effect 
|p9| . For -K^vv the dotted contour giving an ad- 

ditional = 1 is the arc of the circle centered about 
(p, 77) = (1.3,0). It is not quite as restrictive in the SM 
as the 90% CL from the Xg or Atos limit, which is shown 
as the dotted quarter-circle about (p,?7) = (1.0,0). For 
the rest of the SM analysis we drop these three new K 
meson experiments. 

D. Standard Model: (p, r?) Plot 

For the SM (p, 77) plot in Fig. ||, the joint x^ enclosed 
contours are at Icr, 90% CL, and 2cr. The half-circles 
about (p, ?/) = (0, 0) are the center and Icr contours for 
|V^tj*b^d/V^cb^cd|. The hyperbolas are the center and Icr 
contours for e. The quarter circles about (p, r?) = (1,0) 
are for \Vtd\ from Xd in Bd — Bd mixing. The lines em- 
anating from (p, 77) — (1,0) are the central and la lim- 
its for sin (2/3). The Am^ 90% circular arc is shown as 
the dashed quarter-circle, although the analysis weights 
each Ams or each Vtd in X^- We see the effects of the 
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SM Preser; 1 a, 90% CL, 2a 




-O.a -0.6 -0,4 -0.2 0.2 0.4 0.6 0.8 1 



p-T] Plot 

FIG. 5: The (p, 77) plot for the standard model, showing the 
la, 90% CL, and 2a contours of the joint fit, and the central 
and la contours of the various constraints. 

Xs = {Airis/Ts) = 1.35xd\Vts/Vtd\'^ lower bound in the 
SM limiting the length of Vtd oc — p)^ + rf and cut- 
ting off p for /9 < 0. 

E. Four Down Quark Model: {p, rf) Plots 

As in the SM, the plotted p and 77 are taken as the 
coordinates of V*^^, scaling the base of the b — d unitarity 
quadrangle to unity 

p+zi^ = v:,Vudi\v:,v,d\. (30) 

The unitaritity quadrangle is given by 

iCbKd + y,*fcKd + v^^ytd + y^iy^d = o (31) 

where the last term has limits \Ubd/V*f^Vcd\ < 0.15, as 
will be shown later. The near half circles in j ~ S = 613 
{V*i^ = sise*''") at present are due to S14 or 624 (which 
are related) becoming some of the source of the observed 
CP violation in e, so that S13 is less constrained. Then, 
(5i3 can be closer to zero or 180° so that 77 can also be 
small or zero. For projected sin(2a) = +1, 0, or —1, we 
see regions extended beyond the SM regions, which also 
allow 77 to be small. Examining the effect of each new 
K experiment separately, we find that the tt^vv 
result eliminates the large la negative 77 rings from the 
previous analysis Q. 

F. Fraction of the New FCNC Amplitude in e 

In order to display how the FCNC exchange with 
the new phases in Uds can account for the CP violation 
in eK, we plot the ratio of the FCNC contribution to the 
experimental result. In Fig. ^ (eFCNc/kexptI) is shown 
against the phase of V*^^, which is ^13. In Fig. |7[ while 
epcNC cannot account for the entire e result, it can ac- 
count for 60% of it at a la confidence level. 



PDQM, 1 a, 90% CL, 2a 




-1 -0.5 G O.fa i -1 -0.5 0.5 1 

(a) Present P (b) sin(2a) = - 1 P 




(c) sin(2a) = P (d)s]n(2a) = 1 P 



FIG. 6: The (p, 77) plots for the four down quark model from: 
(a) present data, and for projected sin (2q) values of —1, 0, 
and 1. Contours are at la, 90% CL, and 2a. 



PDQM - Present, la, 90% CL, 2a 

1.5 r 

1 .25 P 




50 100 150 200 250 300 350 

6,3 

(EFCNc/lfiexplO VS. 6,3 

FIG. 7: The ratio (epcNc/lecxpt j) of the contribution of the 
FCNC amplitude to eK as a function of the angle Syj,. 

G. Standard Model: (xs, sin (7)) Plots 

Xs is determined in the SM from 

Xs = l.?,^Xd{\Vts\/\Vtd\f- (32) 

The largest error arises from the uncertainty in \Vtd\^ 
which follows from the present 15% uncertainty in 
^/B^fB = 230 ± 35 MeV from lattice calculations Q. 
In the SM, the B factory measurements construct a rigid 
triangle from the knowledge of a and /3, and removes this 
uncertainty in 7 and Xs in the future. 

From present data for the SM (a;s,sin(7)) plot in 
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SM Preser; 1 a, 90% CL, ?c 




0-5 - 

0.4 - 

0.5 - 

0.2 - 

0.' - 



10 2C .30 40 50 00 70 80 

sin(7) vs. 

FIG. 8: The [xs, sin 7) plot for the standard model with 
present limits with contours at la, 90%CL, and 2(t. 



FDQM, 1 a, 90% CL, 2a 




20 40 60 80 100 20 40 60 80 100 

(a) Present (b) sin(2c() = -1^= 




20 40 60 80 100 20 40 60 80 100 



(c) sin(2a) = (d) sin(2a) = 1 



Fig. |, the limits at 2a are 0.56 < sin (7) < 0.99, and 
16 < < 48. Because of the approximately linear re- 
lation between Xs and sin(7), an exact Xs measurement 
(cx l/|Vi(ip) can strongly constrain sin (7) to ±0.07 in 
the SM. 



H. Four Down Quark Model: (ccj,, sin (7)) Plots 

In the FDQM, the sin (7) range goes down to zero at 
1(7, or —0.4 at 2a (Fig. ||), since 77 now goes down to zero 
at la or to —0.2 at 2a where p « —0.5. A larger sin (7) 
range is thus allowed in the FDQM than in the SM. The 
Xs allowed region in the FDQM is 16 to 60 at la or to 
80 at 2(7, which is also larger than the 2a Xg range of 48 
in the standard model. In the FDQM, there is not an 
approximately linear relation between sin(7) and X s as 
there is in the SM. Thus an accurate measurement of Xg 
still leaves a very large region of sin(7) available in the 
FDQM. A subsequent sin(7) measurement will be needed 
to distinguish between the two models. 

I. The Decay Asymmetry from Bs Mixing, 
sin (2</)s) 

(ps is the small angle in the b~ s unitarity triangle given 

by 

+ v,ivts + v:,Vus = 0. (33) 

In Wolfenstein terms this is 

{A\^)x 1+1 X {-A\^-A\'^{p+ir])) + {A\^{p+ir])) x A = 0. 

(34) 



FIG. 9: The (ajs, sin (7)) plots for the four down quark model 
from (a) present data, and (b, c, and d) for B factory cases 
for values of sin {2a) as labeled. 



Then, sin((^s) = A\'^-q/A\^ = A^yy. This is small in the 
Standard Model where sin(2(/)s) = 2A^77 = O.IO77, or at 
2a 

0.030 < sin (2(/),) < 0.060. (35) 

In the FDQM, as seen in Fig. |l0|, -0.2 < sin (2(^3) < 
0.065 at 2(7, and down to zero at la. Here the range con- 
tinues down to zero since 77 can go down to zero. Hence, a 
value of sin (2(/)s) less than 0.03 would signify a deviation 
from the SM. 



J. Fourth Side of the Unitarity Quadrangle Utd 

Unitarity of the b — d columns has four terms, which 
may be written as 

K*bKd + v;,Vcd + v;,Vtd - Um = o, (36) 

since -Um = ^1*6^4^- (We use Udb = U^d)- t^ie uni- 
tarity triangle is scaled by jV^cfc^cdl to make a unit base, 
the complex plot of Um is also so scaled. The length of 
the UM/\V*^Vcd\ side, as plotted in Fig. In], is thus less 
than 0.15, compared to the unit base in the (p, 77) plot, 
and prefers possibly a more vertical direction. The ac- 
curacy of angles and sides of the unitarity triangle must 
and should reach this accuracy for a good test of the SM. 
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FDQM - Present, la, 90% CL, 2ct 



0.08 

c 

CI 0.06 
0.04 
0.02 






Asymmetry s\n(2(p^') vs 



FIG. 10: The (xs, sin (20s)) plot for the Bs asymmetry 
sin {2<j)a) in the four down quark model for present data, with 
contours at la, 90%CL, and 2a. 



FDQM - Present, 1 a, 90% ,2a 



> 
>" 

•\ 0.1 
E 0.05 




-0.15 

-0.15 -0.1 -0.05 



0.05 0.1 0.1 



(lm(Uj/IV,XJ) vs. (Re(U,J/lV,,V, 

FIG. 11: The complex plot of Udb scaled to make it the 
fourth side of the unitarity quadrangle in the b ~ d unitarity 
plot. 



VI. SIZE OF THE MIXING MATRIX 
ELEMENTS AND MIXING ANGLES 



Bound From Z° 



bb 



The weak isovector part of Z'^ bb is reduced by 
(1 - |\^4hP) through 

V = -^il -\Vib\^) + I sin' 9w + lpt, 



A' = -^(l-|V4.n + ^Pt, 



(37) 



[{V'')^ + [A'')^] , (38) 



where Cqcd = 3(1.0385), and pt = 0.0094 for rut = 174 
GeV. Present data and theory give 



R 



pcxp _ 
— 

theory 



0.21642 ± 0.00073 (39) 
0.2158 ±0.0003 (40) 



We note that the |V4b| effect is to decrease Rt, while the 
experiment is about la above the theory. To lowest order 
in 1^46 1 the FDQM effect is 



bb 



std+FCNC 



'nstd 
bb 



= (1 -2.291^46^ 



This gives a contribution to of 

x' = (0.82 + 0.68 X 10^1^461^)^. 



(41) 



(42) 



This is used as a constraint on all angle choices in the 
fit. Taking the 90% CL hmit at x' = (l■64)^ gives the 
bound on \V4b\ from Rb alone of |V4h| w |s34| < 0.035. 



B. 3-D Matrix Element Lego Plot 

From the 2a surface in the 3D space of the magni- 
tudes of the matrix elements involved in the FCNC, Fig. 
[l^ , we can see the limits and ranges of two of the ma- 
trix elements. To 5% accuracy, V4d = —5146"^", and its 
magnitude ranges from 0.035 to 0.085 at 2a. To 10% ac- 
curacy, V4b — — S34, and its magnitude ranges up to 0.020 
at 2a. The third FCNC matrix element, |V4s|, is bounded 
by 0.0004. This requires a fine cancellation between its 



two components in (—5246**^* — si2Si4e*''"), such that 
S24 ~ S12S14 and 624 = Si4 + TT to get the cancelling minus 
sign. This means that there is effectively only one new 
phase, which we may consider as 614. From the cancella- 
tion, S24 ranges from 0.009 to 0.017 at 2a. The cancella- 
tion is to about 1/20 of the value of S24. The third term 
in V4S, 334(523 + 5125136*''^^), then contributes < 0.0009, 
which is the same order as the partly canceling terms. 
The cancellation does not mean fine tuning since one 
could have parametrized V4S by a single angle instead. 
However, the incredibly small size of \V4s\ < could 
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be considered a fine tuning itself. In comparison to the 
SM CKM matrix we should note that keeping the leading 
terms in the real and imaginary parts, Vcs = 1 + iA^X^rj, 
Vcd = -A - iA^X^rj, and Vts = -AX^ - iAX'^-q. So even 
in the standard model there are matrix elements whose 
imaginary parts are as small as 0(A'*), 0(A'''), and O(A^). 

In Wolfenstein terms, \VAd\ ~ S14 « A^, S24 ~ S12S14 ~ 
A^, \VAb\ ~ S34 < AV2, but \Vis\ < A^ The sequence 
may violate the heirarchical expectation from the 3x3 
CKM matrix. 

In the double FCNC Z° exchange amplitude in Bd—Bd 
mixing, via 

Udb = -V:^Vib « Si4e-'^"S34, (43) 

it is only the (5i4 phase in {Udbf' ~ e~2j(5i4 \t^q^ ^g^jj g^jj^^j 
to the SM box diagram term with its phase of (V^^)^- 



C. Phases 

The cancellation in 14^ to make it small requires 624, ~ 
6ii + n. Thus we can display the phases in a two dimen- 
sional plot of Si4 vs. (5i3, as in Fig. |l^. When S13 is in its 
SM range of 40° (sin(7) = 0.64) to 70° (sin(7) = 0.94) 
the SM terms can be dominant and the small FCNC am- 
plitudes allow each Sn equally. For certain values of Su, 
near 80° and 270°, the new physics amplitudes can be 
dominant and S13 can be large, leading to the enlarged 
(p, ry) contours that can reach 77 w and extend beyond 
to Si3 < 200° at 2a. 



D. FCNC Phase Structure 

Using the V4S cancellation structure with S24 = S12S14 
and 624 — Sii + TT, we can rewrite the V^i matrix elements 
in terms of just one phase in the leading terms 

Vm - -sue'''' « si4e*^^^ (44) 
Vib ~ -S34 (45) 

Vis ~ (S24 - Si2Sii)e''^' + S34S23- (46) 

To leading order, it is clear that only the two new phases 
(of which only one is effectively independent) are included 
in the Vii, and therefore in the Uij and in the FCNC 
amplitudes. The SM phase S13 does not appear in the 
leading terms of Uij. 

The FCNC couphngs, using the cancellation relations, 
are 

Uds = -Sl4 [(S24 - S12S14) + S34S23e~*''"] , (47) 
Usb = S34 [(S24 - Si2Sl4)e"'''" + S34S23] , (48) 
Udb = 5145346-**". (49) 

We note that while S14 and S24 are non-zero, the cancel- 
lation in Vis and the ability of S34 to vanish still allow 
all Uij to vanish. 

E. Variable Determination 

In general for the three complex matrix elements Va, 
one would expect three magnitudes and three phases. In 
determining these from the —Uij = V^jVij however, one 
overall phase would not appear experimentally, due to the 
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V*V structure of the Uij. So we can at best determine 
three magnitudes and two phases from the Uij. This 
agrees with the three new angles and two new phases 
introduced in the 4x4 unitary matrix where an extra 
phase has been removed for the definition of the new 
D'^ down quark. Whereas the three Uij may seem to 
contain three real and imaginary parts to be determined, 
they are not independent, since there is one restriction 
between them, namely that the product 



UdsUsbUbd = \VM\^\V4s\^\Vib\'' 



(50) 



is real. So again, we are left with three magnitudes and 
two phases that can be determined by experiments in- 
volving the FCNC amplitudes, which allows us to deter- 
mine the three new angles and two new phases, just from 
low energy experiments involving the Uij . 

With sufficient energy to produce one D quark, the 
angles sa can each be determined separately by the com- 
bined weak production of uD, cD or tD pairs, or from 
the similar decays of the D quarks. 

The cancelation in has related S24 — Asi4 and 624 — 
614 + IT. Thus there are only effectively two independent 
new angles and one new phase to be determined from 
the five independent components of the Uij, leading to 
an overconstrained system. Finding a consistent solution 
is then a test of the FDQM. Of course, if more variables 
are found to be needed, the mixings to five or six down 
quarks would have to be considered. The present fits 
have found non-zero values for S14 and its related S24. 
Yet S34 may still be small or vanish, and the one new 
independent phase is still to be determined, although its 
determination is coupled to that of the CKM S13 phase. 



F. Unitarity Tests on the CKM Submatrix 



For the s column, since \V4s\ < 0.40 x 10 ^, the devi- 
ation from unitarity of the CKM submatrix is bounded 
by 



V^4,P < 0.16 X 10"^ = O.GA^ 



(54) 



For the b column or t row, since |V4h| = S34 < 0.020, 
the deviation from unitarity of the CKM submatrix is 
bounded by 



\V4b\^ < 0.00040 = 0.17A'' « 



(55) 



For the c row, since S24 = S12S14 = Asi4, the deviation 
of the CKM from unitarity is a multiple of the u row 
result from S14 



O.SA^ < |s24p < 3A^ 



(56) 



Finally, the deviation of IV44P from unity is an overall 
measure of mixing to the fourth down quark 



1 - iVi^l' 



'34- 



(57) 



(58) 



The right hand side is dominated by 5^4 giving 
|y44|2 = l-(0.5^3)A^ 

G. Unitarity Quadrangle Completion 

1. b — d Quadrange 



The orthogonality relation between the b and d 
columns of the 4x4 mixing matrix is 



U2b = 0. 



(59) 



Contained in the 4x4 analysis are tests of the unitar- 
ity of the 3x3 CKM submatrix contained in the 4x4 
FDQM mixing matrix. The FCNC couplings Uds, Ust, 
and Udb measure the deviations from orthogonality of the 
columns of the CKM submatrix, in d — s, s — b, and d^b 
projections, respectively. Their sizes will be discussed in 
subsection G under unitarity quadrangles. 

Bounds on the size of the \Vu\, i = 1, 2, 3, bounds 
the deviation from unity of the sum of the squares of the 
three CKM elements in each column 



(51) 



Similarly, for the rows, the |sj4p measure the deviation 
from unity for the sum of the squares of the CKM row 
elements. 

For the d column or u row, since \V4d\ = si4 ~ 0.035 
to 0.085, unitarity of the CKM three elements of the d 
column or u row is off by 



0.0012 < \Vm\^ < 0.0072, or 



O.SA"* < Wa 



4rf 



<3A^ 



(52) 
(53) 



The fourth side of the b — d unitarity quadrangle, scaled 
to make the base of unit length is U^ij/\V*^Vcb\- From 
Fig. 1 1 . we see that the length of the FCNC quadrangle 
side is < 0.15 in the vertical or imaginary direction, and 
< 0.06 in the horizontal or real direction. The sides of the 
b — d unitarity quadrangle can be written in a modified 
Wolfenstcin form as 



^ub^ud 

V^bVcd 

v:t,Vid 



^^db = -Sl4S34e 



AX* {(j) + 111;), and 



V;,Vtd = A\^ [1 - p - 17^ - Xicj) + li,)] 



(60) 
(61) 
(62) 
(63) 
(64) 



'^db with a co- 



where we have introduced (j) + iip into 
efficient to make the scaled quadrangle A independent. 
An example of the scaled b — d quadrangle is shown in 
Fig. |lj. We see that unitarity requires that the length of 
the FCNC coupling side — C/^^ has to be cancelled by an- 
other triangle side to close the triangle, and that occurs 
in Vf*fjVtd having an addition to the SM formula. The area 
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-U^j^X((|)+n|;) 




FIG. 14: The b — d unitarity quadrangle scaled by AX^, with 
sides given as above. 



of the b — d unitarity quadrangle is computed by adding 
the areas of three sub-triangles and a rectangle 

Area(6 - d) = A^X^ [?y + (1 - p)i^] /2. (65) 

We note that if either or both rj and ip are non-zero, 
CP is violated, and the quadrangle has a non-zero area, 
analogous to the SM unitarity triangle result. However, 
as we will see below, the area of the b — d quadrangle is 
different from those of the other unitarity quadrangles by 
the ip term above. 



2. s — b Quadrangle 

The unitarity orthogonality between the s and b 
columns for the s — b quadrangle is 

V:, Vub + Vcb + V;, Vtb -Usb = 0. (66) 



■''°'fDQW\ Present, la, 90%, 2a 

— 0.25 I 



> 




-0.25 -0.2 -0,15 -0.1-0.05 0.05 0.1 0,15 0.2 0.25 

X10-' 

Re(U,,)/lV,,V,J 
(lm(U,J/IV,,V,J) vs. (Re(U,0/IV,,V,j) 

FIG. 15: Contours for the complex FCNC coupling Usb 
scaled by | V^s Vcb|, which is the fourth side of the s — b unitarity 
quadrangle. 



The first term is AX"^ [p — irf) , the second term is AX^ , 
and the third term is —AX^, to leading order. If we scale 
the base to unit length by dividing by AX^ , then the first 
term side is of order 0.02 in length. From Fig. [ij, the 
fourth side of scaled Usb is of order 0.0001, or 0.5% of 
the small third side of the triangle. The enclosed angle is 
then the same as in the SM, (fis = X^i], and the triangle's 
or quadrangle's area is A^X^r]/2. 



3. d — s Qua' 



The orthogonality relation between the d and s 
columns is 



X 10 
,0.25 



ZD 



FDQM - Present, 1 a, 90%, 2a 




KdK. + v:,V,s + V.lVts - Uds = 0. 



(67) 



The largest sides of the d — s unitarity quadrangle are 
of length A, being the first and second terms, and the 
third term is A^ X^ {1 - p + irj) = 0.0004(1 - p -I- i??). The 
fourth side is the FCNC coupling Usd, which is bounded 
in magnitude by 2.5 x 10~^ = A^, as seen in Fig. |l^. Thus 
the FCNC fourth side is at most 6% of the small third 
side. The angle subtended by the small third side is then 
essentially the same as that by the third and fourth sides, 
being (j)^ = A^X'^i]. The triangle's or quadrangle's area is 
also A'^X^T]/2. 



-0.25 -0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 0.25^ 

X 10"' 



ImCUj VS. RefU^J Re(U3 



FIG. 16: Contours for the complex FCNC coupling Uds 
which is the fourth side of the d — s unitarity quadrangle. 
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H. The Sum Rule for the CP violating B Decay 
Asymmetry Angles 

It was shown before Q that as long as the Penguin 
diagrams in the B decays can be neglected, that the sum 
of the CP violating decay angles, even with new physics 
contributions, is tt modulo tt. This can be seen from 
Eqs. @ and where in the sum of (2a + 2/3), the 
Bd mixing phase cancels out in general, regardless of its 
source, and from Eqs. ( |l3| ) and (|T6|), where in the sum 
of (27 + 40s), the phase from Bg mixing cancels out. 
The other tree amplitude decay phases in these equations 
either cancel or sum to the phase of a product of mixing 
matrix elements which becomes a product of absolute 
values squared, with zero phase. This leads to the CP 
violating B decay angle sum rule |^ 

a + /? + 7 + 20s = TT, mod tt. (68) 

VII. CONCLUSIONS FOR ISO-SINGLET DOWN 
QUARK MODELS 

With much new data, it is still the case that FCNCs 
can contribute significantly to Bd — Bd mixing and to 
Bs — Bs mixing, and give contributions with new phases. 
In the FDQM, all sin (2a) are allowed. In the (p, rf) plane, 
the FDQM allows large regions for p < as opposed to 
the p > regions in the SM, and in particular, those 
where rj goes to zero, both with the present data and 
with the projected sin (2a) values. In new physics mod- 
els then, the SM phase 5, or 77, can be smaller, with 
the other phases causing much of the presently observed 
CP violation. In the (a::s,sin(7)) plots in the FDQM, 



all of sin(7) > is allowed at present in contrast to 
sin(7) > 0.55 in the SM, and with no approximately lin- 
ear relation as in the SM. This will require combining 
experimental results of Xg and sin (7), to find out if the 
results correlate to the narrow linear region of the SM 
analysis. The present range for Xg is from 16 to 48 at 
2-(T in the SM, and from 16 to 80 at 2-cr in the four down 
quark model. The b ~ d unitarity triangle, scaled to unit 
base length, has to be measured to an accuracy of 0.15 
or better to begin to limit a fourth side and to verify the 
SM against the FDQM. 

Each Eg generation also contains an iso-singlet or 
sterile neutrino, which may provide a connection between 
the quark and lepton searches for new physics in terms 
of establishing new particle representations. 

The mass of the lightest singlet down quark in Eg could 
be roughly related to the mixing angle by 

6*^4 ~ TOfc/rTiD, and with 1 1^6 1 6*34 < 0.02 (69) 

from combined fits, that gives 

mD> 2500 X mb = 11 TeV. (70) 

Using the single Rb 90% CL hmit of I6I34I < 0.035, which 
is not as strong as the combined fits, gives m/j > 4 TeV. 
The previous analysis |^] gave a lower limit of 1.2 TeV. 
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